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We calculate correction to the critical temperature of a dirty superconductor, which results from 
the local variations of the gap function near impurity sites. This correction is of the order of T c /Ef 
and becomes important for short-coherence length superconductors. It generally reduces a pair- 
breaking effect. In s-wave superconductors small amounts of nonmagnetic impurities can increase 
| the transition temperature. 
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The effect of disorder is of considerable significance for the physics of superconductivity. The well known picture 
for conventional s-wave superconductors includes a strong pair-breaking role of magnetic impurities |jj versus a weak 
effect from time- reversal symmetric perturbations such as nonmagnetic impurities and other lattice defects 0, in 
agreement with the Anderson's theorem H . Generalization to the case of anisotropic Cooper pairing shows that even 
nonmagnetic scatters are pair breakers for a momentum dependent gap. Such a behavior was observed for heavy- 
fcrmion and copper-oxide superconductors and the ci-wave superconducting gap of the latter is now well established 
^ f§- The original Abrikosov-Gor'kov (AG) theory was extended to calculate changes of the transition temperature T c , 
superfluid density, and residual density of states for arbitrary strength of potential scatters in d-wave superconductors 
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On the experimental side, while some measurements on Zn-substituted Y-123 samples show perfect agreement 
with the AG theory j^] , others observe a substantially slower-than-predicted reduction of T c with increasing impurity 
concentration in Y-123 and La-214 systems [Q or with growing ion damage on irradiated samples [||. Various 
theoretical explanations were proposed for this discrepancy including strong-coupling effects [pj or an effect from the 
van Hove singularity p0| . An alternative hypothesis was suggested by Franz and co-workers pd] , who presented 
numerical evidence for the effect of spatial variations of the gap neglected in the AG theory on thermodynamic 
properties for a model two-dimensional d-wave superconductor. First discussions of such effects in conventional 
superconductors with magnetic impurities date to late 60's ]l2] 0. However, they lack consistency with each other 
• and neglect an important effect of the Friedel oscillations. In the present paper we calculate analytically corrections 
to the AG result for the transition temperature resulting from local variations of superconducting order parameter 
s ! ' induced by nonmagnetic and magnetic impurities for various types of pairing. 

Our main results can be summarized as follows. Spatial variations of A, which appear only at very short distances 
of the order of hp 1 from an impurity site, reduce the pair-breaking effect on T c . The correction to the AG result is 
of the order of T c /Ep and becomes important for short-coherence length superconductors such as high-T c cuprates. 
■ In s-wave superconductors small amounts of nonmagnetic impurities can increase the transition temperature, if other 
effects, such as an energy dependent density of states, are small. The effect of spatial inhomogeneity is enhanced in 
I ■ two-dimensional superconductors compared to three-dimensional systems. 

In our treatment of impurities in superconductors we go beyond the standard assumption of a spatially constant 
order parameter QJ^]. The starting point is the linearized gap equation |l^] 

A(k,q) = -T nk,kOG^k' + iq,k+iq)G^(-k' + iq,-k+iq)A(k,q) (1) 

• , u;,k' .k,q 

for a singlet superconducting order parameter A(k, q) (we suppress spin indices) with k and q related to the relative 
and the center of mass coordinates, respectively. The finite-temperature normal state Green's function G"(k, k') 
in an external impurity potential depends on the two momenta. The retarded electron-electron interaction is taken 
in the form V(k, k') = — g/(k)/*(k'), i.e., we assume pairing according to a nondegenerate one-dimensional irre- 
ducible representation of the crystal point group and employ a cut-off lo c , T c < uj c < Ep, in sums over Matsub- 
ara frequencies. Consequently, the superconducting gap in the mixed momentum-space representation is factorized 
A(k, R) = /(k)A(R) and Eq. (|l]) is rewritten for a function of space coordinate A(x) as 

A(x)= 3 T c ^ J dy[^(|x-y|)+if lw (x,y)]A(y). (2) 
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Here, both the gap and the nonuniform part of the kernel Ki u depend on impurity positions Rj (j = 1, . . . , Ni). The 
order parameter is written in the form 

A(x) = (A) 1 ~_^ ) , (3) 

where (A) = y J dxA(x), (F) — y JcfacF(x), and the function F(x) describes relative variations of the gap in space. 
Integrating the gap equation (||) over the volume and assuming the dilute limit ((F) <C 1) we have 

1 = irE fdxdy{K 0u (\x-y\)+K lu ,(x,y)[l-F(y) + (F)}}- (4) 

The terms without F, when averaged over impurity positions, give the standard AG result, whereas the remaining 
members describe corrections due to spatial variations of the gap. In the low impurity concentration limit, the net 
space variation can be considered as a sum of the variations due to each impurity: ifi^(x, y) = Xw=i ^i^(x, y; Rj), 

F(y) = Y^fii ^(y—^j), rii = Ni/V. Substituting these back into Eq. (Q) we notice that all terms Ki u (x, y; Rj)[F(y— 
Rfe) — (^(Rfe))] with j ^ k vanish after averaging over impurity positions. The impurity averaged gap equation has 
the following form 

ln^ + ^)-^ + p) = ^E J dyKU*,y)F(y) , (5) 

where the difference of the di-gamma functions ?/Vs represents the AG result: p — 1/4ttT c t for nonmagnetic impurities 
in an unconventional superconductor and p = 1/2ttT c t s for magnetic impurities in an s-wave superconductor, r and 
t s being potential and exchange scattering times, respectively; N is the Fermi level density of states. In the dilute 
limit the correction to the AG result is expressed via single impurity characteristics only. Similar equations have been 
earlier obtained in Jl2| . Our consideration departs, however, now from previous works. 

In order to express spatial variation of the gap function for a single impurity we return back to Eq. . The 
relation (|^) is now replaced by A(x) = A[l — F(pc)], since (F) vanishes in the thermodynamic limit for one impurity. 
Neglecting variations of A one the right hand side of Eq. (J2j) we immediately arrive at 

F(x) = - 5 T C ]T / dyK luJ (^y) ■ (6) 

This is the so called non-self-consistent approximation for the spatial variation of the gap and its validity for 

short distances from the impurity is justified later. Using symmetry of the kernel we finally get from Eq. (|^) our main 
result: 

T c0 -T c = STag ( n i ) - ^ / dxF(x) 2 . (7) 



gN 

In contrast to previous studies the above equation takes into full account the effect of spatial variations of the gap, 
including the rapid Friedel oscillations of wave-length kp 1 in the vicinity of the impurity. Though Eq. (Q) contains 
only linear in n, corrections, it holds for both magnetic and nonmagnetic impurities and for an arbitrary type of 
pairing. As is seen from Eq. (]?]), the qualitative role of the correction is always to reduce the pair-breaking effect. 
In the case of nonmagnetic impurities in an s-wave superconductor, when 8TxG( n i) = 0, the transition temperature 
increases as a result of spatial variations. 

Nonmagnetic impurities in s-wave superconductors. In order to estimate these effects quantitatively we now turn 
to the calculation of the relative gap suppression -F'(x) in particular cases, considering first a (5-function impurity 
potential U(x) — US(x) in an s-wave superconductor. The 2x2 Nambu Green's function in superconducting state is 
expressed via a t-matrix as G w (x,x') = G^(x, x') + G^(x, 0)CG5(0, x ')> where G5(x,x') satisfies Gor'kov's equation 
without impurity potential, but with corresponding variations of the gap: 

[iw - Vx)r 3 - An] G£(x,x') = *(x - x') - J F(x)Ar 1 G^(x,x') . (8) 

The function F(~x.) has to be found from the gap equation 

F(x) = || Tr [ T i^(x, x)] , 8G U = G„ - G° . (9) 
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Here G° (k) = [iu> — SkT^ — Ati] -1 , Ek is the quasiparticle dispersion, and t$ are Pauli matrices. The expression for 
the i-matrix is t u — Ut 3 [1 — Gf^(0, 0)C/r 3 ] _1 , which we approximate by t u ~ a.T 3 /(irN ) for a = nN U <C 1. 

The short-distance behavior of F(x) can be found approximating G^ by G° • This is exactly the non-self-consistent 
approximation used above. The real space representation for G° is 

G°(x) = -^{(i W + A n )[/ w+ (x)-/ w _(x)] +ifiT 3 [/ w+ (x)+/ w _(x)]} , (10) 
where f u ±(x) = —jrN /(k F x) e ±^ F x e -nx/v F ^ n = ^2 + Substituting it back into Eq. (g) wc find 

- , x sin2/c F a; v-^ e -2^/ t)F 

F x = agN F ttT^ . 11 

A similar result was obtained for a finite-radius impurity by Fetter flqj . Note that the method used in [|l2||l4| l for a 
magnetic impurity would give zero instead of (|iT| ) and, consequently, no effect on the transition temperature. In the 
limit T — > T c one can neglect A and find the following analytic asymptotes 

~ , , sin2fcpa; , ~ ~ irT r 

F (x )=a F for fc^ 1 — 12 

and Fo(x) = 2agNosm(2kFx)/(kpx) 2 e~ x ^ for x > where £0 = vf/(2ttT c ) is the coherence length. The latter 
long-distance asymptote is incorrect since no such rapid oscillations are present in the Ginzburg-Landau regime. The 
correct 1/x behavior is recovered in the self-consistent analysis. However, for the purpose of calculating the correction 
(0) we need to know only the short distance expression (|l^). The function Fo(x) falls of to a negligible value within 
a few Fermi wave-lengths kp 1 . Such rapid relaxation of the gap has also been seen in numerical simulations |]l7| , |l8|| . 
Accordingly, the integral in Eq. (Q) quickly converges and does not depend on the upper limit. Integral's value is 
An 2 a 2 /k F and, hence, the transition temperature of an s-wave superconductor with nonmagnetic impurities is 

Tc = T c0 + -^\n^ , (13) 
° 2r E F A V ' 

where Ao is the zero-temperature gap of a pure system and l/2r = a 2 ni/irNo is the normal state scattering rate. The 
result (|13| ) does not contradict Anderson's theorem, which is valid only for a spatially homogeneous superconducting 
state ||. 

In order to prove the validity of the non-self-consistent approximation we return back to Eq. (|^) and write its 
solution to the first order in F(x): 

G5(x,x') = G«(x-x')-|dyG^(x- y )F(y)Ar 1 G0(y-x') . (14) 

Substituting this into the gap equation and keeping only terms of first order in impurity potential we find after Fourier 
transformation 

F(p){l + | 5 T^;Tr[T 1 G°(p + q)T 1 G°(q)]}=F (p). (15) 

LJ,q 

The difference with the non-self-consistent solution -Fo(p) appears because of the term in the curl brackets with the 
sum over momenta given near T c by 

I. = -V , " 2 + £q 2 £p+q 2 ^ -E^— ( 16 ) 

The short distance behavior of F(x) corresponds to the large-p asymptote for F(p). In this case we neglected 
Matsubara frequencies u in the above equation and found the remaining integral to be I u = m 2 /Ap. The self- 
consistent correction in Eq. ([H]) is gT^^I^ = jgNo/(£op) and is small for £oP > 1- In the corresponding space 
region x < £0 the non-self-consistent solution (|l^) is a good approximation for the gap variations. Numerical data 
show that this may be an accurate approximation even for a strong impurity potential p7[ . 
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Magnetic impurities. We consider a classical spin, which interacts with conduction electrons via an exchange 
Hamiltonian V = J<5(x)5 i <r l , where a 1 are spin Pauli matrices. In order to express appropriately corresponding 
effects, one should include spin indices in Eq. ([!]). Invariance of the singlet superconducting order parameter with 
respect to spin rotations |l9|] suggests, however, that only rotationally-invariant part of the impurity t-matrix produces 
a non-zero effect. Therefore, we can average i u = V[l — G u (0, O)^] -1 with respect to the direction of spin, though 
it does not imply rapid fluctuations of S. The averaged t-matrix is now diagonal in spin indices and is approximated 
as tu = a 2 /(nN ) 2 G^(0,0) for a = JSnN < 1 Q. Assuming a particle-hole symmetry we have G°(0, 0) = 
—TrNo(iuj + Ari)/f2. Following exactly the same procedure as described above we find for short distances x < £ o 



2 2 — cos 2k F x 



F °M= a {kF x Y • (17) 

For pair-breaking magnetic impurities there is an additional non-oscillatory suppression of the superconducting order 
parameter close to the impurity. Note that Eq. ([l?]) [as well as Eq. (|l^)] fails at very short distances x <C kp 1 , where 
an exact band structure becomes important (see, e.g., (l8)). The correction (fjj) can be nevertheless evaluated using 
the expression ( |l7| ) with the lower limit cut-off ~ kp 1 . The result is 



^ = T c0 --^l-5.W^ln^j, (18) 

where l/2r s = a 2 ni/irN^. The qualitative difference in a way how the gap function relaxes near magnetic and 
nonmagnetic impurities has no effect on the integral in (Q), which is always of the order of 1/kp due to rapid 
relaxation. The difference appears, however, in the dependence of the transition temperature correction on the 
strength of perturbation a, showing an additional smallness of the effect for magnetic impurities. 

Nonmagnetic impurities in unconventional superconductors. The real space formalism used above becomes inap- 
plicable in this case, since local variations of the gap function depend on the exact nature of the equilibrium gap 
A(k) = A/(k) and are generally highly anisotropic. Our goal is to show that the order of magnitude of the correction 
to the AG result remains unchanged compared to the case of nonmagnetic impurities in conventional superconductors. 
For that we rewrite all equations in the momentum representation, including <JG w (k, k') = G" (k)£ w G w (k') and the 
gap equation (0) as A(q) = — \gT^ u k /*(k)Tr[riG a) (k + |q, k— ^q)]. Near T c , the non-self-consistent expression 
for the relative variation of the gap is 

agT^ £k+|p/( k - |p)+ £ k-i P /(k+|p) 



ttN ^ j w {^ 2 +e 2 ){uj 2 + e 2 t ) 

We further neglect p in the arguments of /(k± |p) and average |/(k)| 2 over the Fermi surface; this procedure changes 
the numerical factor but not the order of magnitude. After that the integrand becomes the same as for the Fourier 
transform of (|l2|). Therefore, the correction in Eq. (0) will be of the same order of magnitude and finally we have 

T c = T c0 -f (l-f^) (19) 
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with /~3-5; in the BCS theory T c $/Ep w 0.12 a/£ - We now compare this result with the numerical data of Franz 
et al. |pH| . F° r Co /» = 4.7 these authors found a 14% correction to the AG result, which corresponds to / = 5.5. For 
£o/a = 2.5 their result is a 50% correction to the slope of T c {ni) at n, = or / = 10. Experiment on Zn-substituted 
Y-123 superconductor (£o/a ~ 5) found a 40% discrepancy with the AG theory 0], implying / > 10. Somewhat larger 
values for the factor / found experimentally and numerically are related to the two-dimensional (2D) nature of the 
copper oxide materials (see below). 

Two-dimensional superconductors. We consider the case of an s-wave gap. The effect of nonmagnetic impurities 
can be calculated as in the 3D case with the only change that the Green's functions have to be taken in the 2D form: 
/w±(x) = — TriNoHg 1 ^ [±kFx(l ± ft/ Ef) 1 ^ 2 ]- Using the asymptotic form of the Hankel function of the first kind at 
k F x > 1 we have / w ±(x) = N (2ir/k F x y./* e ±i(kpx-z*/A) e -nx/vF ^ Substituting into Eq. @ we obtain for x < £ 

~ . , 2 cos2fcF£ 

F (x ) = -a ? (20) 
7T k F x 

The characteristic feature of the reduced dimensionality is that the order parameter relaxes more slowly to its bulk 
value. Because of this the integral in Eq. (j?j) formally diverges for the above F(x) and becomes dependent on the 
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upper limit. This yields an additional logarithmic factor \n(Ep / Aq) ~ 2-4 to the correction in Eq. ( p~3[ ) enhancing 
effect of spatial variations near impurities on the mean- field T c in 2D. The previous arguments suggest that this result 
must hold for any symmetry of the gap. 

In conclusion, we have calculated a reduction of the pair-breaking effect on the superconducting transition temper- 
ature due to local variations of the order parameter near impurities. Such an effect is important in short-coherence- 
length superconductors and is further enhanced for quasi two-dimensional systems. In singlet superconductors, spatial 
variations of the order parameter near magnetic impurities are less important for T c than near nonmagnetic. Our 
theory explains numerical data for a model d-wave superconductor Q and shows the significance of such effects 
in high-T c cuprates. We also predict an enhancement of the critical temperature of an s-wave superconductor with 
nonmagnetic impurities. This prediction can be tested numerically extending previous model simulations pd|Jl7[ | and 
may have an experimental relevance to Nd2-zCe :r Cu04 and other electron-doped cuprates, which are believed to be 
layered s-wave superconductors. 
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